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Abstract For hydrogen bonded crystals exhibiting proton transfer along hydrogen 
bonds, namely Ol — H- • 02 < — > Ol • • - H — 02, there is a dichotomy of interpre- 
tation consisting in that while the crystal lattice is a quantum object with discrete 
vibrational states, protons are represented by a statistical distribution of classical 
particles with definite positions and momenta at any time. We propose an alter- 
native theoretical framework for decoherence-free macroscopic proton states. The 
translational invariance of the crystal, the adiabatic separation of proton dynamics 
from that of heavy atoms, the nonlocal nature of proton states, and quantum interfer- 
ences, are opposed to statistical distributions and semiclassical dynamics. We review 
neutron scattering studies of the crystal of potassium hydrogen carbonate (KHCO3) 
supporting the existence of macroscopic quantum correlations, from cryogenic to 
room temperatures. In addition, quantum fluctuations calculated for superposition 
states in thermal equilibrium are consistent with measurements of the correlation 
time. There is no temperature induced transition from the quantum to the classical 
regime. The crystal can be therefore represented by a state vector and the dichotomy 
of interpretation must be abandoned. 
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1 Introduction 

The linear formalism of quantum mechanics extrapolated from the level of electrons 
and atoms to that of everyday life leads to conclusions totally alien to commonsense, 
such as Schrodinger's Cat in a superposition of "alive-dead" states and nonlocal ob- 
servables lfT2l . Such conflicts lead to a dichotomy of interpretation consisting in that, 
while at the microscopic level a quantum superposition indicates a lack of definite- 
ness of outcome, at the macroscopic level a similar superposition can be interpreted 
as simply a measure of the probability of one outcome or the other, one of which 
is definitely realized for each measurement of the ensemble Il6l [30l |3T] [32l [33l [341 . 
For open systems, this can be legitimated by decoherence |50J stipulating that an 
initial superposition state should lose its ability to exhibit quantum interferences via 
interaction with the environment. However, since the quantum theory does not pre- 
dict any definite dividing line between quantal and classical regimes, macroscopic 
quantum behavior is possible for systems decoupled from, or very weakly coupled 
to, the surroundings 0. In principle, there is no upper limit in size, complexity, and 
temperature, beyond which such systems should be doomed to classicality. 

For example, it is a matter of fact that defect-free crystals are macroscopic quan- 
tum systems with discrete phonon states at any temperature below melting or de- 
composition. This is an unavoidable consequence of the translational invariance of 
the lattice. However, the dichotomy of interpretation arises for the so-called "pro- 
ton disorder" in crystals containing O-H- • O hydrogen bonds. The coexistence of 
two configurations at thermal equilibrium, say Ol — H • • • 02 and Ol ■ ■ ■ H — 02, has 
been thoroughly investigated in many systems l40l . Although the light mass of pro- 
tons suggests that dynamics should be quantum in nature, semiclassical approaches 
are widely used to rationalize correlation times measured with solid-state NMR and 
quasi-elastic neutron scattering (QENS). Semiclassical protons are thought of as 
dimensionless particles, with definite positions and momenta, moving in a double- 
wells coupled to an incoherent thermal bath. These protons undergo uncorrected 
jumps over the barrier and "incoherent tunneling" through the barrier. In fact, these 
models describe a liquid-like surroundings, at variance with the spatial periodicity 
of the crystal, and strong interaction with the thermal bath is supposed to lead to 
fast decoherence. By contrast, vibrational spectra provide unquestionable evidences 
that the translational invariance and the quantum nature of lattice dynamics are not 
destroyed by proton transfer. Our purpose is therefore to elaborate a purely quantum 
rationale avoiding any mixture of quantum and classical regimes. 

We shall concentrate on the crystal of potassium hydrogen carbonate (KHCO3) 
composed of centrosymmetric dimers of hydrogen bonded carbonate ions (HCO3 )2 
separated by K + entities. At elevated temperatures, the coexistence of two configu- 
rations for dimers is commonly conceived of as a statistical distribution |5] [TT] [45] 
[46). By contrast, systematic neutron scattering experiments measuring a large range 
of the reciprocal space have revealed the macroscopic quantum behavior of protons, 
from cryogenic to room temperatures, and the theory suggests that this behavior is 
intrinsic to the crystal state (15\ [T8l [T9l l24l . The present contribution is a prelimi- 
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nary attempt to elaborate a consistent presentation of experimental and theoretical 
works currently in progress. 

In Sec. |2] we present the crystal structure and the thermally activated interconver- 
sion of dimers. We emphasize why this crystal is unique to observing macroscopic 
quantum effects. In Sec. [3] we show that the adiabatic separation of proton dynamics 
leads to decoherence-free states. Then, we introduce the theoretical framework for 
macroscopic proton states in Sec.|4]and the double-well for protons in Sec. [5] In Sec. 
[6] the calculated scattering cross-section allows us to interpret neutron scattering ex- 
periments in terms of quantum correlations. In Sec. [7] quantum beats arising from 
the superposition of macroscopic proton states in thermal equilibrium are compared 
with the correlation time determined with QENS. In the conclusion, we emphasize 
that the crystal is a macroscopic quantum object that can be represented by a state 
vector. 



2 The crystal structure of KHCO3 
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Fig. 1 Schematic view of the crystalline structure of KHCO3 at 14 K. Dashed lines through pro- 
tons are guides for the eyes, a = 15.06(2) A, b = 5.570(15) A, c = 3.650(8) A, J3 = 103.97(15)°. 
The ellipsoids represent 50% of the probability density for nuclei. 

The crystal at 14 K is monoclinic, space group P2i ja (C^), with four equivalent 
KHCO3 entities per unit cell (Fig. [TJ. Centrosymmetric dimers (HCO^)2 linked by 
moderately strong OH- • O hydrogen bonds, with lengths J?o— o ~ 2.58 A, are well 
separated by K + ions. All dimers lie practically in (103) planes, hydrogen bonds are 
virtually parallel to each other, and all protons are crystallographically equivalent 
(indistinguishable). This crystal is unique to probing proton dynamics along direc- 
tions x, y, z, parallel to the stretching (vOH), the in-plane bending (SOH), and the 
out-of-plane bending (yOH) vibrational modes, respectively. 

From 14 K to 300 K, there is no structural phase transition. The increase of the 
unit cell dimensions and of the hydrogen bond length are marginal, but the popula- 
tion of proton sites changes significantly (Fig. O. Below w 150 K, all dimers are in 
a unique configuration, say L. At elevated temperatures, protons are progressively 
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Fig. 2 Probability densities for protons in the KHCO3 crystal at various temperatures, determined 
with single-crystal neutron diffraction. The ellipsoids represent 50% of the probability density for 
nuclei. 

transferred along the hydrogen bonds to the less favored sites (configuration R) at 
w 0.6 A from the main position. The center of symmetry is preserved and all proton 
sites remain indistinguishable. There is a general agreement that the population of 
the less favored site (or interconversion degree p) is determined by an asymmetric 
double-well potential along the hydrogen bonds [(5] [TT] [T3] [22], but an in-depth ex- 
amination of proton dynamics is necessary to distinguishing statistical disorder or 
quantum derealization. 



3 The adiabatic separation 

Within the framework of the Born-Oppenheimer approximation, the vibrational 
Hamiltonian can be partitioned as 

= ^H + ^t + ^Hab (1) 

where and Jtflt represent the sublattices of protons (H + ) and heavy atoms, re- 
spectively, while %at couples the subsystems. For OH- • O hydrogen bonds, cou- 
pling terms between OH and O- • O degrees of freedom are rather large 11381 l40l . 
and beyond the framework of the perturbation theory. Two approaches, either semi- 
classical or quantum, are commonly envisaged. 

In the semiclassical view, protons are thought of as dimensionless particles, with 
definite positions and momenta, moving across a potential hypersurface [3] |4] [23] 
[42] 04]. Complex trajectories involving heavy atoms lead to mass renormalization, 
and to incoherent phonon-assisted tunnelling |5] [TT] |4lJ. This approach is quite 
natural when the Born-Oppenheimer surface is calculated from first principles, but 
quantum effects can be severely underestimated. 

Alternatively, if the classical concept of "trajectory", totally alien to quantum 
mechanics, is abandoned, adiabatic separation of the two subsystems, namely M{\ 
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and J^ t , may lead to tractable models HI] [2B [23j [37j |40l |47l . Then, light protons 
in a definite eigen state should remain in the same state in the course of time, while 
heavy atoms oscillate slowly, in an adiabatic hyperpotential depending on the proton 
state, through the coupling term. This separation is relevant for KHCO3 because 
adiabatic potentials for different protons states do not intersect each other. Then, 
protons are bare fermions and quantum correlations should occur ||T8l . 

In fact, the separation is rigorously exact in the ground state, since protons should 
remain in this state for ever, if there is no external perturbation. Furthermore, for 
asymmetric double-wells, with wave functions largely localized in each well (see 
below Sec. [5]i, the adiabatic separation should also hold for the lowest state of the 
upper minimum and long-lived superposition states should interfere. 



4 Macroscopic proton states 

Consider a crystal composed of very large numbers N a , N\,, N c (jV = N a N\,N c ) 
of unit cells labelled j,k,l, along crystal axes (a), (b), (c), respectively. The two 
dimers per unit cell are indexed as j,k,l and j',k,l, respectively, with j = f. For 
centrosymmetric dimers, there is no permanent dipolar interaction, so that inter- 
dimer coupling terms and phonon dispersion are negligible ll22l l25l l26l . The eigen 
states of the sublattice of protons can be therefore represented in a rather simple 
way with the basis sets of eigen states for isolated dimers 

A H1-H2 dimer is modelled with coupled centrosymmetric collinear oscillators 
in three dimensions, along coordinates aijki and 0C2jki (<* = x,y,z), with respect to 
the center at Oojki- The mass-conserving normal coordinates independent of j,k,l, 
and their conjugated momenta, 

a s = — = (at - a 2 + 2ao), P sa = —j= (Pia - Pia), 

Y Y (2) 

a a = — (ai + a 2 ), Paa = — |(^la+^2a), 

lead to uncoupled oscillators at frequencies h(o sa and h(D aa , respectively, each with 
m = 1 amu. The difference (h(o sa — tXOaa.) depends on the coupling term (say X a ). 
The wave functions, f^j kl (o!a), ^^((Xs — V20q), cannot be factored into wave 
functions for individual particles, so there is no local information available for 
these entangled oscillators. Consequently, the degenerate ground state of indistin- 
guishable fermions must be antisymmetrized. For this purpose, the wave function is 
rewritten as a linear combination of those for permuted oscillators as 



©ojki± = -4 EI %jki(«a) f^ojki^ - ^2«o) ± ^ojki(« s + V2ao) 
V2 a L 



(3) 



and the antisymmetrized state vectors as: 
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\0jkl+)®\S) = |o 0jkl+ )®-L[|ta 2 >-UiT2)]; 



\0jU-)®\T) = |0 O jki-)®-^ 



I TiT 2 > + 1 iil 2 > + ^[! Ti4.2> + 1 iiT 2 >] 



(4) 



The oscillators are now entangled in position, momentum, and spin. In contrast to 
magnetic systems (9), there is no level splitting, so the symmetry -related entan- 
glement is energy-free. It is also independent of X a - Furthermore, contrariwise to 
Keen and Lovesey (29), or Sugimoto et al. |43l . we argue, as an experimental fact, 
that there is no significant exchange integral for protons separated by w 2.2 A ifPTll . 
Neutron diffraction and spectroscopy show that protons in KHCO3 are neither de- 
localized nor itinerant particles and there is no sizable energy band structure. 

In quantum mechanics, normal coordinates (O define nonlocal pseudoprotons 
(m = 1 amu), say <^ s jki an d ^ a jkh with an internal degree of freedom corresponding 
to symmetric or antisymmetric displacements of two "half-protons", respectively. 
Each dimer site is a superposition of two such half-protons. Obviously, pseudo- 
protons are totally alien to the intuitive conception of particles, based on classical 
mechanics, but they are the actual observables, whereas individual particles are not. 

Consider now the sublattice of protons. The spatial periodicity leads to collective 
dynamics and nonlocal observables in three dimensions. With the vibrational wave 
function for the unit cell j, k, I, namely Sqjkix = ©Ojkix ± ©oj'kft' where T = "+" or 
"— ", usual phonon waves can be written as 

j N c N b N a 

S()x ( k ) = ~r^7 EEL s ojki* ex p( ik • L ) ' ( 5 ) 



^ l=lk=lj=l 

where k is the wave vector and L = ja + kb + /c, with the unit cell vectors a, b, 
c. This equation represents collective dynamics of H-H dimers thought of as com- 
posed bosons. This would be correct if the crystal was composed of indistinguish- 
able dimer entities (KHCC»3)2. However, neutron diffraction shows that the crystal 
structure is composed of KHCO3 entities related to each other through the appropri- 
ate symmetry operations. The probability density of each atom is equally distributed 
over all equivalent sites and, conversely, the probability density at each site includes 
contributions from all indistinguishable atoms of the same kind. Consequently, the 
sublattice of protons must be thought of as a sublattice of nonlocal indistinguishable 
fermions and antisymmetrization of the plane waves §5$ leads to 

kL = 0modulo27T. (6) 

Consequently, there is no phonon (no elastic distortion) in the ground state and 
this symmetry-related "super-rigidity" [T8j is totaly independent of proton-proton 
interaction. Then, the lattice state vectors in three dimensions can be written as: 



|So+(k = 0)>®|S>; 
|S _(k = 0)}®|T}. 



(7) 
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Each macroscopic state of the sublattice represents a nonlocal pseudoparticle with a 
mass m = 1 amu, namely a pseudoproton, 3? a or S? s , with a definite spin-symmetry 
and an occupation number of per site. There is no local information avail- 

able for these entangled states and the wave functions So T (k = 0) represent collec- 
tive oscillations of the super-rigid lattice as a whole with respect to the center of 
mass of the crystal. Finally, the ground state of the sublattice is a superposition of 
the pseudoproton states as: 



This ground state is intrinsically steady against decoherence. Irradiation by plane 
waves (photons or neutrons) may single out some excited pseudoprotons. Entangle- 
ment in position and momentum is preserved, while the spin-symmetry and super- 
rigidity are destroyed. However, the spin-symmetry is reset automatically after de- 
cay to the ground state, presumably on the time-scale of proton dynamics. Conse- 
quently, disentanglement reaches a steady regime such that the amount of transitory 
disentangled states is determined by the ratio of the density-of-states for the sur- 
rounding atmosphere and external radiations, on the one hand, and for the crystal, 
on the other. This ratio is so small that disentangled states are too few to be ob- 
served, but they allow the super-rigid sublattice to be at thermal equilibrium with 
the surroundings, despite the lack of internal dynamics. The main source of dis- 
entanglement is actually the thermal population of excited proton states. However, 
even at room temperature, the thermal population of the first excited state (< 1 % for 
yOH w 1000 cm -1 ) is of little impact to effective measurements. 

For the sublattice of bosons in the isomorph crystal of KDCO3, (O and (O are 
not relevant. There is neither spin-symmetry nor super-rigidity. Dynamics are rep- 
resented with normal coordinates (f2]i and phonons (01. Needless to say, the H and 
D atoms have the same number of 12 degrees of freedom per unit cell, but the sym- 
metrization postulate shrinks the size of the allowed Hilbert space from ~ 12*" for 
bosons to ~ YLJV for fermions. 



5 Proton dynamics 

The interconversion degree at thermal equilibrium (Fig. |2j is determined by the po- 
tential function for protons. On the one hand, the bending modes do not play any 
significant role, since they show rather modest anharmonicity and the population of 
excited states is negligible. On the other hand, the double-wells for the OH stretch- 
ing (Pig. [3j is known from experiments. The distance between minima (2x ~ 0.6 A) 
is given by the crystal structure. The upper states at ZiVo2 and ZiVo3 were determined 
from infrared and Raman band profiles |[T3l[T4l . and the ground state splitting (/?Voi) 
was observed with incoherent inelastic neutron scattering (HNS) ll22ll26l . The po- 
tential obtained through best fitting exercises is over determined and largely model 




0))®|S>; 

0)>®|r>. 



(8) 
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L (A) R (A) 

Fig. 3 Potential function and wave functions for the OH stretching mode along the hydrogen bond 
intheKHC0 3 crystal. V(x) = 374x + 0.4389 x 10 6 + 5516exp(-30.8;t 2 ). V andx are in cm" 1 and 
A units, respectively 1 1 3 11221 . The oscillator mass is 1 amu. 

independent. In addition, the oscillator mass of 1 amu is not a free parameter. It is 
determined by 2xq for a given set of energy levels. 

For the |0) and 1 1) states, the potential asymmetry leads to substantial localization 
of the wave functions in the lower and upper wells, respectively. However, tunneling 
is possible through the tiny delocalized fraction (e « 0.05) visible in Fig. [3] 

This potential has been a puzzle ever since it was determined because the up- 
per minimum was naively thought of as corresponding to the transfer of a proton. 
However, this is unlikely, for this would lead to unrealistic dimers composed of 
di-protonated (H2CO3) and non-protonated (CO^) entities l20l . Such entities are 
ruled out by the centrosymmetric character of proton dynamics established by the 
symmetry -related selection-rules observed in the infrared and Raman 1361 . It is now 
clear that, if pseudoprotons are the observables, this nonlocal potential accounts for 
pseudoproton dynamics along x a or x s . The |0) < — ► |1) transition corresponds to 
the through-barrier transfer (tunneling) of a pseudoproton as a rigid entity, with no 
energy transfer to the internal degree of freedom. The IINS bandwidth, very close 
to the spectrometer resolution [22], shows that this transition is virtually dispersion- 
free and there is no visible splitting suggesting any difference for the transfer of 
or £P S l22l l26l . On the other hand, the upper states |2) and |3) correspond to 
excitations of internal stretching coordinates, v a (infrared) or v s (Raman). They are 
slightly different, but this is unimportant for interconversion since thermal popula- 
tions are strictly negligible for these states. 

The interconversion dynamics involving 8P a and &* s can be rationalized with the 
potential surface along coordinates x a and x' s — x s ± v2*r> 




Fig. 4 Schematic view of the tunneling wave functions. For the sake of clarity, the weak compo- 
nent of the wave function in one dimension is multiplied by a factor of 2. 



r(x a ,4) = y(^)+y(x' s ). (9) 

The energy level scheme comprises three states at 0, hvoi (with twofold degeneracy) 
and 2h Voi- The non- symmetrized local wave functions (Fig.|4|i are simple product 
of the local wave functions in one dimension (Fig. [3]). Proton configurations for the 
three states are tentatively sketched in Fig. [5] The ground state corresponds to the 
structure observed at low temperature, with both £P a and ^ s in the L configuration. 
The antisymmetrized macroscopic state analogous to (O, namely |0+)|S}|0— )|7*}, 
can be obtained via (01 to Similarly, the upper state vector at 2/iVoi with both 
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Fig. 5 Schematic view of proton configurations for the tunnelling states. In the ground state LL 
and in the upper state RR at 2hVo\ , protons are fully entangled. In the intermediate state LR (or 
RL) at hvoi all proton sites are equally occupied. The spin symmetry and the super-rigidity are 
destroyed. 

pseudoprotons in the R configuration is |1+)|S}|1— }\T). In the intermediate state 
at IiVqi, only one pseudoproton (either ^ a or is transferred to the R configura- 
tion, so the spin-symmetry and the super-rigidity are destroyed. Then, plane waves 
(O lead to state vectors |l+,0-,ki ) and |0+, 1 — ,k i>. Note that the |0) < — > |1) 
transition is effectively observed with IINS, thanks to energy and momentum trans- 
fer, whereas the |0) < — > |2) transition cannot be probed directly, according to the 
quantum theory of measurements. 

The proton transfer degree calculated supposing the three levels at thermal equi- 
librium is 

p (T) = [ P0l (T) + 2p 2 m (T)} [1 +p i(T) + poi (T) 2 }- 1 , (10) 

where poi(T) = exp(—hVo\/kT) is the probability for the transfer of a pseudopro- 
ton. The dashed line in Fig. [6] clearly shows that this equation is at variance with 
observations. 

In fact, owing to the adiabatic separation of proton dynamics, energy exchange 
with the surroundings occurs exclusively via photons, with the momentum conser- 
vation rule kio + koi = 0. Consequently, re-entanglement occurs spontaneously in 
the intermediate state as 

|i+,o- ; k 10 ) + |o+,i- ; koi)=2- 1 /2[|o + )|5)|o-)|r) + |i+)|5)|i-)|r)] (ii) 

and the transfer degree 



p(T)=2p 2 01 (T)[l+p 2 0l (T)}-\ 



(12) 
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Fig. 6 Temperature effect 
on the proton transfer degree 
in the KHCO3 crystal. Stars: 
experimental according to 
fill . Solid line: eq. flfj for 
two-levels. Dashed line: eq. 
HO) for three-levels. 




50 100 150 200 250 300 350 
T(K) 



is in reasonably good agreement with measurements (see the solid line in Fig. 
The energy difference of 2/jVoi between RR and LL configurations is therefore con- 
firmed and Figure [6] is an indirect evidence that the degenerate intermediate state 
observed with IINS is not thermodynamically stable, thanks to the purely quantum 
re-entanglement mechanism ( TPTV 

However, the interpretation of Bragg diffraction is ambiguous because the LL 
and RR configurations of the proton sublattice are crystallographically equivalent, 
as they are related through a translation vector (a/2,b/2,0) (see Pig . |5j . The recip- 
rocal lattices are identical and it is unknown whether neutrons were diffracted by 
either sublattice, with probability 1 — p and p, respectively, or by a superposition 
state, ( 1 - p ) 1 1 1 10+) \S) |0-) | T) + p '/ 2 1 1 +) \S) 1 1 -) | T) . The former case is a mix- 
ture of LL and RR configurations analogous to disorder under consideration in many 
crystallographic |27l|45j|46), solid-state NMR J5][39l and QENS HQ works. Alter- 
natively, a superposition should give rise to quantum interferences corresponding to 
coherent fluctuations of the probability density at proton sites. 



6 Probing quantum entanglement with neutrons 

Neutrons (spin 1/2) are unique to observing the spin-symmetry of macroscopic 
states @. However, quantum entanglement is extremely fragile, because it is not 
stabilized by any energy. Consequently, only "noninvasive" experiments, free of 
measurement-induced decoherence, are appropriate |34| . For neutron scattering, 
this means (i) no energy transfer (ii) no spin-flip and (iii) particular values of the 
neutron momentum transfer vector Q preserving the super-rigidity. (By definition, 
Q = k, — k/, where k, and k/ are the initial and final wave vectors, respectively.) 

The dotted lines in Fig. Q] enhance the network of double-lines of proton sites in 
dimer planes. We present below neutron scattering experiments for (i) double-lines 
of protons, (ii) arrays of double-lines in two dimensions, and (iii) the sublattice in 
three dimensions. For the sake of simplicity, it should be born in mind that elastic 
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scattering events are identical for LL and RR configurations and, therefore, indepen- 
dent of the interconversion degree. 



6.1 Double-lines of entangled protons 



Consider an incoherent elastic neutron scattering (IENS) experiment conducted with 
(i) the best resolution in energy, in order to reject inelastic scattering events, (ii) a 
modest resolution in Q, so Bragg peaks merge into a continuum and long-range cor- 
relations are overlooked. For momentum transfer Q a along a, the scattering func- 
tion for an entangled pair can be written as: 



II 

| (OjklT, | expiga («2jki - aqjki) + TjT f expiGa (« 2 jki + «o jk i) \0jkln) 
x (OjklTf\ expiga («ijki - «ojki) + T f Ti expiQa (aiju + a jki) \0jkln) \ 2 
Xexp(-2Wi a )5(ffla)- 



(13) 



Each bracket represents a scattering event by a pseudoproton located at both sites 
(iaojki). The product of two brackets means that each neutron is scattered si- 
multaneously by the two pseudoprotons superposed at the same sites, either in- 
phase, |±) < — > |±) (TfT; = +1), or anti-phase, |±) < — > (TfT; = —1). The spin- 
symmetry is probed along the neutron-spin direction with 100% probability and 
there is no spin-flip because the initial and final states are \0jklTi) and \0jklt{), 
respectively, for one scattering event and vice-versa for the other one. Thanks to 
adiabatic separation, the lattice Debye-Waller factor exp(— 2W^a) can be factored. 
The energy transfer is hco a and 8(co a ) accounts for energy conservation. In the har- 
monic approximation, the scattering function is Ifl31l24l 



S±± (Ga,«*x) 
COS 4 (2 a Oo) 



exp 



Qa u 0a 



S±^(Qa,COa) = 
sin 4 (g a ao)exp 



Vl+4A a 



exp-e a M 0a 



-Go 



*0a 



'0a 



exp(-2WLa)5(C0a), 
exp(-2W La )5(co„). 



(14) 



Here, Uq u = %j (2m(Oo a ) is the mean square amplitude for uncoupled harmonic os- 
cillators in the ground state. The intensity is proportional to the incoherent nuclear 
cross-section for protons, Oni ~ 80.26 b (1 barn = 10~ 24 cm 2 ), and the gaussian-like 
profiles for uncorrelated scatterers are modulated by cos 4 (Gee c«o) and sin 4 (Gee c«o) 
(Fig. 0a). 

Such interference fringes were effectively observed with the MARI spectrom- 
eter UJ at the ISIS pulsed-neutron source (Fig. [7] b) |24]. Best fit exercises yield 
Oo-values in reasonable accordance with the crystal structure and the estimated os- 
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cillator mass is virtually equal to 1 amu. These experiments are positive evidences 
of pseudoproton states with spin-symmetry (0). 




Fig. 7 (a): Comparison of the theoretical profiles S(Q, 0) for a non-entangled pair (dot-dashed), 
and for an entangled pair of fermions, according to lilt (solid line). The doted and dashed curves 
represent interference fringes for in-phase and out-of-phase scattering, (b): S(Q y ,0) measured at 
20 K for a single-crystal of KHCO3 (solid curve with error bars). Comparison with the best fit 
(dashed line) obtained with I14t convoluted with a triangular resolution function. The dashed line 
with error bars is the difference spectrum. *: Triangular functions due to other scattering events. 



In fact, neutron plane waves are scattered coherently by double-lines of entangled 
pairs perpendicular to Q a , which are reminiscent of Young's double-slits. However, 
the interference fringes ([Pil l are clearly different from those anticipated for dis- 
tinguishable (classical) double-slits, which should be proportional to cos 2 (g a ao) 
iflOl |48l |49l . Equation (TBI is also at variance with the scattering function for 
protons delocalized in a symmetric double-wells. In this case, there is no definite 
spin-symmetry for the tunneling states |0+) (ground state) and |0— ) (at TicOt). The 
scattering functions for elastic scattering |0±) < — > |0±) and inelastic scattering 
|0±) < — > |0=f) are 

S±±(Qa,(O a ) = cos 2 (e a ao)exp(-e2 M 2 a _2W La ) d(co a ), _ 
S±^ (Qa,(Oa) = sin 2 (2a«o)exp (-Qiul a -2W La ) 8{co a - ©t). 

Then, interferences evidence that a single proton is located in two wells. They are 
visible if the instrument can effectively resolve the tunnel splitting. Otherwise, com- 
plementary fringes would merge into the gaussian profile anticipated for a single- 
well. 

Clearly, scattering by a superposition of entangled double-lines with spin cor- 
relations cannot be confused with other double-slits experiments. The fringes are 
evidences of nonlocal pseudoprotons and there is no means whatever to probe the 
local particle behavior. 

Since ( fT3l l and ( TBI hold in the same way for LL and RR configurations, inter- 
ferences are independent of the interconversion degree. However, at elevated tem- 
peratures, the intensity at large Qa-values is depressed by the lattice Debye-Waller 
factor and fringes are less visible. 
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Needless to say, interferences are neither expected, nor observed, for KDCO3 
CD- 



6.2 Diffraction 

A necessary condition for noninvasive neutron diffraction is that Q x , Q y , Q z , should 
match a node of the reciprocal sublattice of protons, so neutrons probe super-rigid 
states without any induced distortion. The only information conveyed by such events 
is the perfect periodicity of the sublattice, so the Debye- Waller factor is equal to 
unity at any temperature. In addition, thanks to the spin-symmetry, the scattered 
intensity is proportional to the total cross-section Oh ~ 82.0 b lfl8l[35ll . Otherwise, 
the spin-symmetry is destroyed, so the intensity scattered by protons is proportional 
to the coherent cross-section Ohc ~ 1-76 b and to the Debye-Waller factor for non- 
rigid lattices. The enhancement factor (Th/Ohc ~ 45 is quite favorable to observing 
quantum correlations. Furthermore, the intensity scattered by the sublattice of heavy 
atoms, proportional to (7cKC0 3 ~ 27.7 b, is depressed by the Debye-Waller factor 
exp— 2Wkco 3 (Q)- Therefore, the contribution of heavy atoms at large Q-values is 
rather weak, compared to that of the entangled sublattice, especially at elevated 
temperatures. 

The dashed lines in Fig. [TJshow that proton sites are aligned along x and y, but 
not along z- Consequently, the noninvasive condition can be realized for Q x and 
Q y exclusively, whereas Q z never coincides with a nod of the reciprocal lattice of 
protons. Then, the diffraction pattern depends on whether we consider incoherent or 
coherent scattering along Q z . 



6.2.1 Super-rigid arrays in two dimensions 

Consider diffraction by super-rigid arrays in (103) planes and incoherent scattering 
along Q z . In the unit cell, there are two indistinguishable double-lines parallel to y, 
so the periodicity of the grating-like structure is D x /2, withD x wa/ cos 42° w 20.39 
A. On the other hand, the spatial periodicity of double-lines parallel to x is D y = b. 
The differential cross-section for a superposition of pseudoproton states (O can be 
then written as 



hrM 



K N b 

EE! [ ex Pi2y (kDy-yo) + TjTfexpiQy (kD y +y )] 



j=ik=i 2 (16) 

x [expigx (JDx/2-xq) +TiT f expi(2 x (JD x /2 +x )}} 2 exp-2W z (£ z ), 

with N' a = Neutrons are scattered either in-phase (TfT; = +1) or anti-phase 
(TfTj = — 1) by orthogonal pairs of lines separated by 2xo as 0.6 A and 2yo ~ 2.209 
A, respectively. The phase matching condition, namely xq (yo) commensurable 



Macroscopic Quantum Entanglement 



15 



- ■ 



k^O * . ■ 




m 

[5 -10 -5 


5 SO ] 



■15 -Ifl -5 5 lfl IS -15 -IH 



^3 



•f; » 



lA') 



k=6 


« 
- 

. s : 

... r n 

■ *J 




* -5; 


*■ 







-15 -10 -SO 3 10 i* -H -HI 4 




k = 2 

-: : ='i 4 i! 


1p k 













k - 5. 





5 10 -5 1 

lA 




k = S 

■It. 

i 




^ .IS.' 
V 




I -10 -s 


S 10 1 

"'> 


11 





ID 35 -15 -10 -5 5 H 15 -15 -lfl -5 5 10 15 

(A" 1 ) (A 1 ) 



Fig. 8 Cuts of the diffraction pattern of KHCO3 at 300 K in various (a*,c*) planes. The arrows 
emphasize ridges of intensity parallel to Q z and perpendicular to the dimer plane (dash lines along 

GO- 



with D x /2 (D y ), is intrinsic to the crystal structure. {•••} 2 accounts for simul- 
taneous scattering by the superposed pseudoproton states, without neutron-spin 
flip. The compound Debye-Waller factor exp—2W z (Q z ), including contributions 
from all atoms, corresponds to incoherent scattering along Q z . The permutation 
xq < — ► —XQ,yo < — ► —yo gives the same equation for RR and LL configurations. 
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Table 1 Orders n y , n x , and positions Q y , Q x , of rods of intensity arising from the entangled 
array of orthogonal doubles lines of protons in two dimensions. Obs. Q x : positions in A -1 of the 
observed rods of intensity along Q z in Figs[8]and[9] Q y D y /n is rounded to integers. 



"v 


Qy (A"') 


QyDy/n 


Tiff 


k = Q y /b* 


Gx 


Obs.fi* (A- 1 ) 











+ 1 





n x 7t/xo 


0,±10 


1 


2.86 


5 


-1 


2.57 


(n x + \/2)n/xo 


±5,±15 


2 


5.71 


10 


+ 1 


5.14 


n x n/xa 


0,±10 


3 


8.57 


15 


-1 


7.71 


(n x + \/2)n/x 


±5,±15 



The diffraction pattern is composed of rods of diffuse scattering parallel to Q z , 
cigar-like shaped by the Debye-Waller factor, at Q x , (2 y -values corresponding to 
divergences of (fl6l) , Such divergences occur at Q y = n y 7t/yo ~ « y x 2.86 A -1 (Table 
[TJ, since Q y D y /n w 5« y is integer. Contrariwise, there is no divergence for anti- 
phase scattering at Q y = ±(rc y + 1 /2)n/yo, because Q y D y /n « 5(« y + 1 /2) is not 
integer. For n y even, Q y D y / % is also even, Tj = Tf, and ridges are anticipated at 
gx = n x n/xo k n x x 10 A , since Q K D x /n f=a 68« x is even. Alternatively, for n y 
odd, Q y D y / % is also odd, Tj 7^ Tf, and ridges are anticipated at Q x = (n x + 1 /2)n/xo, 
since QxD x / n « 68 (n x + 1 /2) is even. 




The cigar-like shaped rods were effectively observed with the SXD ||2][28] instru- 
ment at the ISIS pulsed neutron source (Figs[8]-[T0]l. For k = 0, they appear at Q x = 
and ±(10.00 ±0.25) A -1 , in accordance with 2xq k 0.6 A. For k = 1, they are barely 
visible. For k = 2 or 3, we observe ridges at Q x = ±(5 ±0.2) and ±(15 ±0.2) A" 1 , 
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still along Q z . These features are best visible in Fig. [9] for k = 2.6, in accordance 
with Table Q] There is no visible ridge at k = 4. Then, from k = 5 to 9, we observe 
in Fig.[8]the same sequence as for k = to 4, and rods at k = 7.7 in Fig. [9] 



Fig. 10 Diffraction pattern of 
KHC0 3 at 30 Kin the (a*, c*) 
reciprocal plane at k = 0. The 
arrows emphasize ridges of 
intensity parallel to Q z and 
perpendicular to the dimer 
plane (dash lines along Q x ). 




Similar diffuse scattering was observed at low temperature in the k = plane, 
at 14 K Ull and 30 K (H) (Fig. [lOj. As anticipated, the rods are unaffected by 
the interconversion degree. However, at low temperatures, they are partially hidden 
by the anisotropic diffuse intensity, centered at Q = 0, due to elastic and inelastic 
incoherent scattering by protons. This continuum precludes observation of the ridges 
at ±5 A -1 for 2 < k < 3. Quite paradoxically, quantum correlations are best visible 
at elevated temperatures. 

As anticipated from Sec. HI the same experiments performed with a crystal of 
KDCO3 do not evidence any cigar-like shaped ridge of enhanced intensity, in addi- 
tion to regular Bragg peaks, for the sublattice of bosons |fl9l . 



6.2.2 Super-rigid arrays in three dimensions 

From ( fT6l ). the differential cross-section for the LL or RR pseudoproton states is 
written as 



d£2 .q tf 



N' a N b N c 

ILL {[expigyftoy-.vo) +TiT f expig y (^£) y +yo)] 

j=l k=l 1=1 



[expigx C/At/2 -x ) +T;Tfexpigx (jD x /2 +xo)]expiQ z lD z } 2 



2 



(17) 



This equation describes no spin-flip scattering events that do not destroy the spin- 
symmetry. Divergences occur along the previous rods of intensity at Q z = ±n z 2n / D z , 
with D z Ri c x cos28° w 3.28 A and 2n/D z m 1.92 A -1 . These enhanced peaks are 
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visible in Figs[8lfT0l even at rather large (2 z -values, thanks to super-rigidity. At el- 
evated temperatures, they clearly emerge from the rods of diffuse scattering ( [ToT l 
depressed by the Debye-Waller factor. These enhanced peaks were not observed for 
the deuterated crystal |[T9l . 

Experiments presented above, in this present section, are clearly consistent with 
pseudoprotons forming decoherence-free macroscopic single-particle states with 
spin-symmetry. The underlying theoretical framework presented in Sec. 2]is based 
on fundamental laws of quantum mechanics. There is no ad hoc hypothesis or pa- 
rameter. The adiabatic separation, clearly validated by observations, can be regarded 
as an intrinsic property of hydrogen bonds in this crystal. 



7 Quantum interferences 

Superposition of decoherence-free proton states must lead to quantum interferences, 
or quantum beats. The non-antisymmetrized wave functions for the states |0) and 1 1) 
(Fig- E]l can be written as ||20l : 

fojki = cos0 yo(x-x m ) +sin0 yb(*+*m); 
fijki = - sin0 \l/o(x-x m ) +cos0 \j/ (x+x m ); 

where x stands for x a or x / s and y/Q(x±x m ) are harmonic eigen functions for 
the second-order expansion of the potential around the minima at ±x m ; tan20 = 
v 0t/ ( v 0i — v 0t)> where /iVot « 18 cm -1 is the tunnel splitting for the symmetric po- 
tential. Then, cos as 1 and sin (j) = e sw 5 x 10~ 2 . Superposition leads to harmonic 
oscillations of the probability density at the beating frequency Vob = 8eVoi ~ 4vot ~ 
2.5 x 10 12 , in proton per second units (H s" 1 ) lf8l [T6l l20l . 

The LL < — > RR fluctuation rate can be rationalized with two distinct mecha- 
nisms, either single-step or two-stepwise (Fig. 0J. The single-step mechanism cor- 
responds to superposition of the states at and 2/tVoi, corresponding to LL and RR 
configurations (Pig. |5j, respectively. The interconversion rate due to quantum beats 
is: 

Vib = 2evDb«p(-2Av i/*r). (19) 

For the two-stepwise process, firstly, either ^ a or ^ is transferred at k = to 
the LR configuration (Fig. |5j with probability exp(— ItVoi/kT). Secondly, this state 
undergoes fast re-entanglement (fTTT i leading to the upper state (configuration RR) 
with probability exp(— 2liVoi/kT). The interconversion rate is then 

v 2b = 2vbbexp(-3AvDi/AT). (20) 

These theoretical rates must be compared to QENS measurements of a KHCO3 
crystal, from 200 to 400 K ifTTI . The scattering geometry (Q || x) was selected in 
order to probe proton dynamics specifically along the hydrogen bonds. The inverse 
relaxation time (or attempt frequency), t^ 1 = 2 x 10 12 s _1 , is sufficiently close to Vob 
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to suggest that (i) QENS and vibrational spectroscopy techniques probe the same dy- 
namics and, (ii) the two-stepwise mechanism < f20b is prevailing, in accordance with 
the larger pre-factor. In addition, the measured rate follows an Arrhenius law with 
an activation energy £ a = (336 ± 32) cm -1 significantly different from 3/iVoi ~ 648 
cm" 1 . In fact, d2"0t accounts for coherent fluctuations of two pseudoprotons, with 
probability exp(— 3hvoi/kT), and pre-factor 2vot>, while QENS probes the fluctua- 
tion rate of a pseudoproton with probability exp(— 3hvo\/2kT), and pre-factor Vob- 
Hence, 3/zVoi/2 = 324 cm -1 accords with E a . It transpires that the QENS technique 
is an incoherent probe of coherent oscillations of the proton probability, because 
neutrons are plane waves, rather than a particle-like probe of incoherent stochastic 
jumps |fTT|. 

In fact, semiclassical models J5] QT| ED are based on inappropriate premises, (i) 
The potential asymmetry supposedly due to static effects of neighbouring dimers 
should be temperature dependent. This is at variance with the interconversion de- 
gree ( fT2b and Fig. [6] (ii) Coupling to phonons is posited to be necessary to mediate 
the through-barrier proton transfer at low temperatures. This is not relevant within 
the framework of the adiabatic separation, (iii) A smooth transition to the Arrhenius 
behaviour of classical jumps is supposed to occur at elevated temperatures. Con- 
trariwise, neutron diffraction shows that there is no transition to the classical regime 
(Figs[8l-[T0il. In addition, $1% and d2Qb show that an Arrhenius behavior is not nec- 
essarily an evidence of the semiclassical regime. Logically, these incorrect premises 
lead to confusing the activation energy E a with the potential barrier FTP . 



8 Conclusion 

It is often argued that a complex system in continuous interaction with its environ- 
ment should be in a significantly mixed state that cannot be represented by a state 
vector. In marked contrast to this widespread opinion, we have accumulated con- 
sistent experimental evidences that the sublattice of protons can be represented by 
a state vector at any temperature up to 300 K. This macroscopic object exhibits all 
features of quantum mechanics: nonlocality, entanglement, superposition and quan- 
tum interferences. There is no transition to the classical regime because the plane 
waves of the thermal bath cannot destroy entanglement intrinsic to the lattice pe- 
riodicity. The spin-symmetry of proton states can be transitorily destroyed but the 
decoherence degree is insignificant because the density-of-states of the surroundings 
is negligible compared to that of the crystal. 

The cornerstones of the theoretical framework are: (i) adiabatic separation; (ii) 
the fermionic nature of protons; (iii) indistinguishability and degeneracy. There is 
no ad hoc hypothesis or parameter. Entanglement is intrinsic to the crystal symme- 
try, irrespective of the strength of proton-proton interactions. Dynamics is rational- 
ized with pseudoprotons forming macroscopic single-particle states with remark- 
able spin-symmetry and super-rigidity. These quantum correlations are effectively 
probed with neutrons and quantum interferences arising from entangled double- 
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lines or long-range correlations in three dimensions emphasize that protons in the 
crystal field are not individual particles possessing properties on their own right. 
Collective dynamics suggest that the whole crystal should be conceived of as a mat- 
ter field that is a superposition of macroscopic single-pseudoparticle states. In addi- 
tion, super-rigidity adds a new item, along with superfluidity and superconductivity, 
to the list of quantum "super" properties in the condensed matter. 

The interconversion degree at thermal equilibrium is consistent with the double- 
well potentials for pseudoprotons determined from vibrational spectra and quantum 
beats accord with QENS measurements of the fluctuation rate. The double-wells 
for protons is invariant over the whole range of timescales from vOH vibrations 
(~ 10~ 15 s) to diffraction, through QENS, and at any temperature. 

This review emphasizes that the dichotomy of semiclassical protons in a quan- 
tum crystal lattice should be abandoned. There is every reasons to suppose that this 
conclusion holds for many hydrogen bonded crystals and macroscopic quantum be- 
haviors open up new vistas for further investigations. 
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